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ABSTRACT 

We calculate the average stretch-corrected rise-time of type la supernovae (SNe la) in the Supernova 
Legacy Survey. We use the aggregate lightcurves of spectroscopic and photometrically identified 
SNe la to fit the rising part of the lightcurve with a simple quadratic model. We obtain a lightcurve 
shape corrected, i .e. stretch-corrected, fiducial rise-time of 17.02^Q'2|(stat) days. The measured 
rise-time differs from an earlier finding by the SNLS (|Conlev et al.l [2006h due to the use of different 
SN la templates. We compare it to nearby samples using the same methods and find no evolution 
in the early part of the lightcurve of SNe la up to z — 1. We search for variations among different 
populations, particularly subluminous objects, by dividing the sample in stretch. Bright and slow 
decliners (s > 1.0) have consistent stretch-corrected rise-times compared to fainter and faster decliners 
(0.8 < s < 1.0); they are shorter by 0.57lg'gg(stat). Subluminous SNe la (here defined as objects with 

s < 0.8), although less constrained, are also consistent, with a rise- time of 18-03±5 j^stat). We study 
several systematic biases and find that the use of different fiducial templates may affect the average 
rise-time but not the intrinsic differences between populations. Based on our results, we estimate 
that subluminous SNe la are powered by 0.05 — 0.35M Q of 56 Ni synthesized in the explosion. Our 
conclusions are the same for the single-stretch and two-stretch parameterizations of the lightcurve. 



1. INTRODUCTION 

Type la supernovae (SNe la) are believed to be very 
bright explosions of carbon-oxygen white dwarfs (WDs) 
in binary systems that undergo thermonuclear runaway. 
Despite this consensus, the nature of the companion and 
explosion mechanism are under debate. The enormous 
energy release of a SN la takes around 2 weeks to peak in 
the optical. During this early phase of the SN, when the 
luminosity rises to maximum, fundamental aspects of the 
explosion, and possibly of the progenitor, are revealed to 
the observer. The spectroscopic presence of intermediate 

*Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Sci- 
ences de PUnivers of the Centre National de la Recherche Sci- 
entifique (CNRS) of France, and the University of Hawaii. This 
work is based in part on data products at the Canadian Astron- 
omy Data Centre as part of the Canada-France-Hawaii Telescope 
Legacy Survey, a collaborative project of NRC and CNRS. 
jgonzalez@astro.utoronto.ca 



mass elements (between carbon-oxygen) such as silicon, 
besides iron-peak elements, suggests incomplete burning 
from a subsonic deflagration compone nt in the explosion 
(|Filippenkolll997HMazzali et alj|2007t) . A subsequent de- 
layed detonation is typi cally favored to account for the 
energy output of SNe la (|Khokhlo"v| [T991; Yamaok a et al.l 
119921 (Wooslev fc Weaverlll994D . 

The radioactive decay of the 56 Ni synt hesized in the 
explosion powers the lightcurv e of SNe la (jTruran et al.l 
119671 IColgate fc McKed 119691) . It is thought that the 
variety in the amount of 56 Ni gives rise to most of the 
scatter in SN la luminosities. Values of 0.4 — O.8Af 
of 56 Ni result in "normal" SNe la. "Subluminous" or 
1991bg-like objects, very faint and fast evolving SNe la, 
are produced by small amounts of 56 Ni: 0.07 — 0.1 7M Q 
(jFilippenko et alj|1992t iLeibundgut et al.lll993t ). Slowly 
evolving and very bright SNe, s uch as 199 IT are g ener- 
ated by 0.85 - l.LM of 56 Ni (jPhillips et al.l [19921) . In 
this framework, the duration from explosion until maxi- 
mum light, known as the rise-time, is partly determined 
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by the amount of 56 Ni fe.glArnet t 1982; Shigc vama et al 
l99 l | Brancbj|1992t iKhokhlov et al.l 1 19931: iHoeffich et al 



1993; IPoim'nguez et al.ll200lD .' 

The picture of the early behavior of SNe la has 
greatly benefited from the advent of improved aspheri- 
cal and m ulti-dimensional simulations of delayed detona- 
tions (e.glBravo fc Garci'a-Senzll2008l: IKasen et al]|2009l : 
iMaeda et al.ll2010t iRoepke et al.l l2010f K These authors 
show that the quantity of 56 Ni synthesized in the explo- 
sion that powers the lightcurve depends primarily on a 
single parameter that can be characterized by the den- 
sity at which the supposed transition from deflagration 
to detonation in the delayed detonation framework oc- 
curs. If the amount of burning in the deflagration phase 
is very large (due to, for example, numerous symmet- 
ric sparks), the transition density is low and the deto- 
nation produces less 56 Ni. Nevertheless, other factors 
like the meta llicity could affect the luminosity to sec- 
ond order (e glTimmes et al.ll2q03HGallagher et al.ll2008t 
Howell et all 120091: iNeill et al.l 120111: iBravo et al.ll2010D . 



Hoflich et al.l (|2010l ). for example, predict that a depen- 



dence of the early lightcurve on the C/O ratio of the 
progenitor (determined mainly by the progenitor main- 
sequence mass) potentially affects the slope of the rise 
behavior at epochs earlier than 5 days before maximum. 

Beyond the physics of the explosion, the rise-time 
study of SNe la can directl y prov ide information on the 
progenitor system. IKasen! (|2010f) suggests that in the 
single degenerate (SD) scenario, where a WD accretes 
material from a main-sequence or a red giant, the early 
phase of the lightcurve can be affected by the interac- 
tion of the companion star with the SN shock, particu- 
larly at shorter wavelengths. The magnitude depends on 
the observing angle, as well as on the size and separa- 
tion of the companion, and can be directly tested with 
well sampled early rise-time data. In t he double detona - 
tion of sub-Chandrasekhar mass WDs, iFink et al.l (l2010h 
predict faster rise-times for brighter and slowly declining 
SNe la. They also show that models with additional ra- 
dioactive decay of other short-lived nuclei like 52 Fe and 
48 Cr contribute and magnify the early epoch output of 
SNe la. 

These models suggest that the lightcurve parameteri- 
zation of SNe la could need additional corrections besides 
the standard width-luminosity -and color -luminosity as 
well as the recent host mass-luminosity ([Sullivan et al.l 
l2010HKellv et al.ll2010fl . especially at early epochs. The 
inclusion of secondary parameters that correlate with lu- 
minosity may provide an improved calibration of SNe la 
as distance indicators for cosmological studies. On the 
other hand, secondary parameters might not affect the 
luminosity directly but still be very important to under- 
stand the SNe la themselves. 

Subluminous SNe la are an extremely faint population 
with redder colors. Standard scenarios do not reach the 
very faint explosions of subluminous SNe la and cannot 
reprod uce their Fe- and Si-composition. iPakmor et al.l 
(2010) propose a merger system of two WDs of similar 
0.9M,7,) that is able to account for the low 56 Ni 
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mass of ~ O.IMq. These objects are useful to understand 
the amount of 5 °Ni needed to power a lightcurve, and to 
investigate possible secondary effects. As subluminous 
SNe la are preferentially found in older, more metal- 



rich stellar populations (lHamuv et aT1ll996t lHowellll2001l : 
IGonzalez-Gaitan et al.ll201 If) . thev also offer an opportu- 
nity to study metallicity effects on the amount of 56 Ni 
produced. 

The early epoch study requires good early sampling of 
SN lightcurves. Modern rolling surveys like the Sloan 
Digital Sky Survey (SDSS) and the Supernova Legacy 
Survey (SNLS) have a good cadence and provide very 
well sampled pre-maximum lightcurve data. Compar- 
ing fiducial stretch-corrected rise-times measured from 
different groups is an extremely difficult task, partic- 
ularly because of the rather arbitrary definition of a 
fiducial, or, s = 1 template. Ideally, the study of 
several SN la populations, be it at varying redshift or 
host environment for instance, should be done with the 
same techniques and template s. Previou s studi es in- 
clude the pioneering work by iRiess et al.l (|T999) who 
found a B-band fiducial rise-time of 19.5 ± 0.2 days 
using combined near b y da ta —10 days before maxi- 
mum. lAldering et al.l ([20001 ) find no evide n ce of evolu- 
tion in SN la rise-time and IConlev et al.l (|2006l here- 
after C06) confirmed this by calculating the rise-time of 
the low-z data compared with higher- z from the SNLS 
(z < 0.9) finding consistent rise-times, 19.58+°'^ for the 
nearby and 19.101° ^(sta tistical)±0.2(systematic) for 
the SNLS. IStrovinkl ([2007ft found a fiducial rise-time of 
17.44 ± 0.39 days with an algorithm that avoids external 
templates and fits smooth functions to the photometric 
data directly. Furthermore, they found individual rise- 
times in a bimodal dist ribution aroun d 18.81 ±0.36 days 
and 16.64 ± 0.21 days. lHavden et all (|2010b[ ) propose a 
new parameterization of the lightcurve: they divide it 
into a rising and a falling component, and find a stretch- 
corrected rise-time of 17.38 ± 0.17 da ys for the SPSS , 
rejecting the bimodal distribution of IStrovinkl ([2007) . 
They also find evidence for a faster stre t ch-corr ected 
rise-time in slow d ecliners. lHavden et al.l (]2010aD and 
IBianco et al.l (1201 ID go further to constrain the SD mod- 
sn (2010) and have largely ruled out red giant 
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companions. 

None of the previous studies have included sublumi- 
nous SNe la. In this work, we make use of the large, 
homogeneously observed, but intrinsically heterogeneous 
SN la sample obtained by the SNLS, in particular of 
the phot ometric group of low-st r etch (s < 0.8) objects 
found in IGonzalez-Gaitan et al.l ([201 ID . We investigate 
the stretch-corrected rise-time behavior of different SN la 
stretch populations to look for any secondary variations 
in the lightcurve besides the standard width-luminosity 
relation. We study the importance of different templates 
and lightcurve fitting techniques on the stretch-corrected 
final rise-time calculations. For proper comparison with 
other samples, we also analyze nearby and SDSS data 
with our own techniques. The paper is organized as fol- 
lows: In § [2] we present the rise-time parameterization 
used in our study, in §[3] we introduce the data we used, 
its analysis follows in § 2] and comparisons between sub- 
luminous SNe la and other populations, as well as to 
previous studies are presented in § [5] We conclude in 

2. RISE- TIME FIT 
2.1. Lightcurve parameterization 
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We follow a similar approach to C06 to measure 
the common stretch-corrected rise-time of our combined 
SN la sample data. We call this stretch-corrected rise- 
time a "fiducial" rise-time. Different SNe la are cor- 
rected for lightcurve (LC) shape and normalized to peak 
flux, after expressing all data in the same filter - nor- 
mally B - via if-corrections. For the ^-corrections, we 
use the appropriate spectral energy distribution (SED) 
based on the epoch, SN's redshift and LC parameters 
such as stretch and color. For restframe epochs (cor- 
rected by dividing by (1 + z)) before explosion, where the 
flux is generally zero, the i^-corrections on ly include the 
redsh ift term correction: — 2.51og(l + z) ( Nuge nt" et al.l 
2002). Due to the small number of low-stretch SNe la 
available, our study of the rise-time focuses on the av- 
erage properties of the aggregate stretch-corrected LCs, 
instead of individual rise-times. 

The LC fitter used throughout th is study is a modi- 
fied version of t he SiFTO algorithm (jConlev et al.ll2008t 
IGuv et al.f2010D that uses SED templates to model SN la 
LCs. The SED sequence has been trained on a la rge sam- 
ple of low-z and SNLS data (IConlev et al.ll2008fl and de - 
scends from the spectral sequence of Ifisiao et al.l (|2007l) . 
Although the heterogeneous low- z training sample does 
not include subluminous nor SN2002cx-like objects, the 
modified fitter also works we ll for subluminous objects 
(|Gonzalez-Gaitan et alll2011l hereafter Gff). The LC 
is parameterized with a single time-axis stretch param- 
eter in the restframe B-band, sb- Each observed filter 
is stretched differently according to a relat ion between 
the wavelength and sb ()Conlev et al.ll2008h . The same 
SEDs are used to if-correct all data to the restframe. 
Observed epochs are transformed to effective dates, r, 
by r = it — t())/s(l + z), where to is the date of maxi- 
mum flux in the restframe -B-band and z is the redshift. 

The LC fits normally use all available photometry, in- 
cluding the early time. In order not to affect the rise- 
time analysis, we fit the LC and calculate stretch with 
only photometry past a certain phase, i.e. r t < r < 35 
days. We call this the "core" LC. The transition value 
from rise-time to core LC, r*, considered here ranges from 
— 10 to —8 days (see § 14. 5p . One can also use other LC 
parameter s based uniquely on post-maximum data, such 
as Am is (iPhillipsI 119931: (Phillips et al.l ll~999l Alterna- 
tively, lHavden et al.l (|2010bl hereafter H10) suggest a 
two-stretch approach for dividing the LC fit around the 
epoch of maximum light To. If fit) is the flux, we have: 

f(A = ffo F ((t-to)/s r ) if * < to m 

J[l> \f F((t-t )/s f ) if* >*o { ) 

where fo is the flux at peak (which is the same for both 
parts of the LC), Fit) is the fiducial (s = 1) flux of the 
template, and s r , Sf are the rise and fall stretches, respec- 
tively. This means that for the two-stretch technique, the 
epoch is given by r = (t — to) /s r (1 + z) for pre-maximum 
data and by r = (t — io)/s/(l + z) for post-maximum 
data. In this paper, we use both the single- and two- 
stretch parameterizations to fit the LCs of SNe la. 

2.2. Rise-time parameterization 

To study the rise-ti me, we adopt the d escrip - 
tion used by iRiess et all (|1999ft : lAldering et~aT1 (|2000l ): 
iGoldhaber et all (|2001[ ) and C06. If r is the LC phase 



and t r is the fiducial rise-time of a s = 1 SN that we 
want to constraint, we have for the flux: 

f(r) = a{T + t r ) n + b, (2) 

for t > t r and b at earlier epochs. Normally, the nuisance 
parameter b should be zero, but we investigate possible 
photometric subtraction errors that could systematically 
shift fluxes up or down (see § I4.6[) . a is the slope of 
the rise and is another nuisance parameter in the fit. 
The quadratic form, n — 2, is motivated by the charac- 
teristic expansion of a hot object in the Ray leigh- Jeans 
tail of the SN la SED at early times in the B- and V- 
filters. The flux is proportional to the temperature and 
to the radius squared, the latter in a homologous expan- 
sion (r cx vt) depending linearly on the velocity of the 
expansion: / cx r 2 T cx v 2 (t +t r ) 2 T ', where r is the radius, 
v the rise "speed" and T the temperature. Since AT/T 
is almost constant observationally, and Av/v is also close 
to constant at early epochs, then / oc (r + i r ) 2 . Although 
past studies agree that this quadratic model fits the rise- 
time behavior quite well, small variations among SN la 
species are possible, not only for the rise-time t r but also 
for the power n of the rise. We therefore investigate such 
variations as well (see § 14.71) . 

2.3. 56 Ni radioactive decay 

The radiated luminosity of a SN la is powered by the 
radioactive decay of 56 N i through 56 Co to 56 Fe (jArnettl 
119821: lArnett et ail 119851: IBranchl fl992[) . At peak, this 
energy is still trapped in nickel and is slo wly released by 
its radioactive decay ()Nugent et a,l.lll995h : 

Lboi = ~fS(T r )M m , (3) 

where Mni is the mass of 56 Ni, 7 ~ 1.2 is the ratio 
of bolometric to radioactivity luminosities for normal 
SNe la, and S is the luminosity from radioactivity per 
unit nickel mass given by: 

S= (6.31 e - r '~ /8 - 8 + 1.43 e- T "/ m ) x lO^ergs^Afg 1 . 

(4) 

This comes from the characteristic decay times of 56 Ni 
and 56 Co of 8.8 and 111 days, respectively, with mean en- 
ergy release per decay of 1.71MeV and 3.76MeV. There- 
fore, by measuring the individual rise-time r r and bolo- 
metric luminosities L& ; of SNe la, we can estimate the 
amount of 56 Ni synthesized in the explosion under the 
assumed model. Variations in the amount of 56 Ni mass 
produced in the explosion will lea d to changes in the ris e- 
time for a fixed peak brightness (jContardo et al.|[200Q[ ). 

The value of 7 is near unity but could in principle 
vary slightly depending on the 56 Ni mass, i. e. for large 
nickel masses, the radioactivity energy released above 
the photosphere will escape directly instead of remain- 
ing trapped, leading to lower values of 7. In the absence 
of a model for its v ariation, we assume 7 = 1.2 ± 0.2 
([Nugent et al.l fl995T ) and propagate this uncertainty in 
our calculations. 

1 Note that we use t r for the average fiducial rise-time and r r 
for the individual rise-time of each SN uncorrected for stretch 
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3. DATA 
3.1. SNLS 

The five-year (2003-2008) Supernova Legacy Sur- 
vey took advantage of the MegaCam wide- field imager 
(jBoulade et al.l [20031 ) to observe supernovae, as part of 
the Canada-France-Hawaii Telescope (CFHT) Deep Syn- 
optic Survey. It used a rolling search technique suited for 
full coverage of the LC, including early-time, by imaging 
every 3-4 days (2-3 days in SN restframe) four fields in 
four different bands <7mTm»m-Zm , which are si milar to the 
filter system ugriz of the SDSS ([Smith et alJ l2002). The 
SNLS had two independent p hotometric recon structions 
based in Canada and France (jGuv et al.l 120101 ) . We use 
the Frenc h photometry he r e, as in the SNLS co smological 
analyses ([Guv et alj|20lot iConlev et alJ[20TTh . 

C06 used 73 spectroscopically confirmed normal SNe 
la for their rise-time study. Here, we use 30 photometri- 
cally identifie d low-stretch (s < 0.8) SNe la and 221 nor- 
mal SNe la ([Gonzalez-Gaitan et al.l 120 lit iPerrett et al.l 
1201 If) at0.1<z<0.7to study the respective stretch- 
corrected rise-times and compare them. The photomet- 
ric candidates went through a selection process that en- 
sured proper redshift and stretch estimation. As part 
of the selection criteria to constrain the stretch, we re- 
quire here good lightcurve sampling, particularly enough 
early data. For this, we re-fit for stretch using only data 
after r t , to ensure the lightcurve shape (stretch) is not 
affecting our rise-time estimation. We require at least 
two data points within r t < r < days of restframe B- 
band maximum. The restframe £?-band for these objects 
will be mostly mapped by tm (and somewhat by 
The lightcurves must also have a good SiFTO fit with 
a total reduced chi-squared, xtt better than 3.0. Addi- 
tionally, we require that the error in the determination 
of the peak epoch be less than 0.7 days (1.2 days for the 
two-stretch parameterization, as shown in section [4.2jl . 
For the rise-time region fit, SNe la are required to have 
data prior to —10 days. These cuts reduce the sample to 
134 SNe, of which 12 have s < 0.8. 

In our dataset, 35% of the low-s and 74% of the normal 
SNe la have a spectroscopic redshift. The effect of having 
uncertain photometric redshifts can lead to considerable 
errors propagated in the fC-corrections as well as in the 
effective dates r of each SN, ultimately affecting the rise- 
time fit. Throughout this study, photometric redshifts 
are used in the same way spectroscopic redshifts are, but 
a proper inclusion of their errors is analyzed in § 14.11 

3.2. Low-z and SDSS 

To test our techniques, we also analyze available low-z 
rise-t ime data (see Ta b le [B, as well as the SDSS-II SNe la 
from iHoltzman et al.l ((2008) . We can thus compare our 
results to low redshift to study any evolutionary effects. 
We cannot provide a comparison of our low-s stretch- 
corrected rise-time analysis to nearby ones, as published 
early data for nearby s < 0.8 are scarce. The cuts applied 
to the SNLS are also used for the SDSS-II and low-z 
data, i. e. at least two data points in restframe B-band 
at n < t < days, and at least one point before r t . 

4. ANALYSIS 
4.1. Rise-time calculation 



To make a combined rise-time study of all SNe la in 
e ach sample, we process each lightcurve as described in 
§ 12.11 We then take the restframe B-fluxes of the low- 
stretch SNe la and overlay them on a single plot, as 
shown in Figure [TJ For the rise-time region, we use early 
data, before a characteristic transition epoch r t = —10 
days to fit it to the model of Equation [5] (see Figure [5]) . 
Minimizing the \ 2 of the best fit model with respect to 
the data provides the best stretch-corrected rise-time. 
The choice of the transition epoch used affects the fit, as 
will be shown in § 14.51 We perform an iterative outlier re- 
jection (similar to C06) to discard any data points that 
strongly deviate from the fit according to Chauvenet's 
criterion (|Taylorlll997[ ). For the total sample, this corre- 
sponds to 3.4cr and discards 25 (out of 449) data points, 
but does not affect the final fit. 

The early time fit of the aggregate LC, nevertheless, 
does not take into account uncertainties in the SN fit pa- 
rameters and the photometric redshift. Errors in th e red- 
shift estimate induce small changes in stretch (jGllh . Ad- 
ditionally, the use of a common SED for the LC fit as well 
as for the if-corrections introduces strong correlations, 
as do the uncertainties of maximum date and stretch of 
each SN fit. We handle these by shifting the photomet- 
ric redshift according to its error (for those SNe la with 
photo-z) in a Monte Carlo (MC) fashion. For each new 
redshift a corresponding fit provides the new LC parame- 
ters (stretch, date of maximum and flux scales) and their 
respective errors (and covariance matrix). These values 
and errors, in turn, are used to generate a new set of 
random LC values (for SNe with and without photo-z) 
that determine the model SED for each epoch. The K- 
corrections are recalculated for the new set of redshift 
and LC parameters. For each of these realizations we fit 
the rise-time data of all SNe to the model and get the 
best fiducial, stretch-corrected rise-time t r and slope a, 
respectively. The final values are drawn from the median 
of all realizations combined. The statistical uncertainties 
in the parameters are given by the la uncertainties in the 
distribution (see Figure [3]). The distribution is broader 
for low-s because of the smaller sample size, and because 
the low-s sample includes more SNe with photometric 
redshifts, and hence larger errors in z (see Gil for de- 
tails). 

Final rise-times for the SNLS, SDSS-II and nearby 
samples are shown in Table [5J The different samples 
agree well with each other providing no evidence for evo- 
lution. This is confirmed by the fiducial rise-time of 
the SNLS SNe la at z > 0.7. Rise-times for different 
SNLS stretch populations after the primary LC-width, 
or stretch, correction are also consistent. However, the 
data are also in agreement at a ~ la level with a sec- 
ondary slight trend for longer fiducial rise-times in fast 
decliners. High-stretch SNe (s > 1.0) have the short- 
est stretch-corrected rise-times, ~ la shorter than mid- 
stretch SNe (0.8 < s < 1.0), and low-stretch (s < 0.8) 
also seem to follow this trend. In the next sections, we 
investigate systematic effects that could modify these re- 
sults. 

4.2. Comparison of single- and two-stretch fits 

Here, we investigate differences arising from the pa- 
rameterization of the LC in two pieces (rising and falling) 
as in Equation 12. 11 and as done in H10. We process our 



The Rise-Time of Subluminous SNe la 



5 



TABLE 1 

Nearby objects used in the rise-time calculation. 
Name (column 1), redshift (column2), single-stretch (column 3), 2-s rise- and fall-stretch (columns 4 and 5), and source 

(COLUMN 6) ARE PRESENTED. 



Name 


Redshift 


Stretch 


Rise-Stretch 


Fall-Stretch 


Source 


SN1990n 


0.0034 


1.07 


1.18 


0.96 


Lira et al. ( 19981 


SN1994ae 


0.0043 


1.05 


1.09 


1.03 


Riess et al. (2005l,Altavilla et al. (20041 


SN1997bq 


0.0094 


0.91 


1.05 


0.81 


Jha et al. ( 20061 


SN1998aq 


0.0037 


0.97 


1.10 


0.90 


Riess et al. (2005) 


SN2000e 


0.0048 


1.06 


1.28 


1.00 


Valentini et al. (20031; Tsvetkov (2006a1 


SN2001el 


0.0039 


0.96 


1.12 


0.91 


Krisciunas et al. (20031 


SN2001v 


0.0150 


1.09 


1.22 


1.05 


Hicken et al. (20091 


SN2002bo 


0.0042 


0.94 


1.02 


0.93 


Krisciunas et al. (2004) 


SN2003du 


0.0064 


1.00 


1.04 


1.00 


Anupama et al. f2005l,Stanishev et al. (20071, Leonard et al. (20051 


SN2004fu 


0.0092 


0.90 


0.99 


0.84 


Tsvetkov (2006b1; Hicken et al. (20091 


SN2005cf 


0.0064 


0.98 


1.03 


0.96 


Pastorello et al. (2007l.Wansr et al. (20091. Hicken et al. (20091 


SN2005mz 


0.0176 


0.63 


0.91 


0.53 


Hicken et al. (20091 


SN2006cp 


0.0223 


1.05 


1.11 


1.02 


Hicken et al. (20091 


SN2006gz 


0.0237 


1.22 


1.30 


1.19 


Hicken et al. (20071 


SN20061e 


0.0174 


1.08 


1.10 


1.05 


Hicken et al. (20091 


SN2007af 


0.0055 


0.94 


1.07 


0.91 


Hicken et al. (2009) 





Fig. 2. — Overlaid B-band restframe lightcurves of low- (s < 0.8) and normal-stretch (s > 0.8) SNe la from the SNLS in the early rise 
at r t < —10 days. Red squares are outliers from the fit (~ 3.4<r). The solid black line is the best quadratic fit, the blue dotted line is the 
median of the flux points, the dashed violet line is the fiducial Gil template used, and the red dot-dashed is the C06 template. 
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ply due to the increased complexity of the model, or 
if there is a real advantage in using it. To compare 
the two models, we use the Bayesian Information Cri- 



17 18 

t r (days) 

Fig. 3. — Distribution of 100 MC stretch-corrected rise-time fits 
to s < 0.8 (filled orange) and s > 0.8 (empty blue) SNe la. The 
dashed and dotted lines represent the median for both samples. 



data by fitting it with a 2-s model, and we find that the 
rise-stretch and fall-stretch are consistent within the er- 
rors, although the average rise-stretch is slightly larger 
than the fall-stretch: (s r - s f ) = 0.09 (RMS of 0.24). 
The mean rise-stretch is also larger than the stretch 
from the l-s model, (s r - s) = 0.06 (RMS of 0.17) (see 
Figure 2]). This can be explained as the day of maxi- 
mum gets pushed to later times by more than half a day, 
(to(s r , s f ) - t (s)) = 0.65 (RMS of 1.77). 

The mean error in the peak date is much larger for the 
2-s fits, (Ato(s r ,s f ) - A* (s)) =0.87 (RMS of 0.53) (see 
Figure [5]). An incorrect estimate of the date of maximum 
has a direct effect on the rise-time calculation. For our 
MC technique, in which the fit parameters are varied 
according to their errors, large uncertainties in the peak 
epoch will affect the final rise-time even more. Our initial 
cut of At < 0.7 days rejects more than 80% of the 2-s 
fits of the SNLS sample (40% for the l-s fits) , including 
all low-s objects and all z > 0.7 SNe la. We therefore 
relax the cut to Ato < 1.2 days for the 2-s fits. 

We do not change our low-s definition for the 2-s fits 
to s r < 0.8 as this would reflect only the early part of 
the LC. If we were to use the fall-stretch instead, some 
low-s candidates would have even lower stretch values. 
Instead we note that most candidates have a mean 2- 
s stretch, i. e. (s r + s/)/2, closer to the single-stretch: 
(s r + s/)/2 - s = 0.01 with 0.05 RMS, and use this 
to separate the different stretch samples. This definition 
keeps most of the low-s candidates and is the best overall 
description of the lightcurve. 

With the new stretch values we derived for the 2-s fits, 
we repeat the procedure to estimate the fiducial rise- 
time, using the rise-stretch instead of the single stretch 
values for the ^-correction and normalization of all prc- 
maximum epochs. We obtain shorter stretch-corrected 
rise-times than the l-s approach, about 1-2% for the dif- 
ferent SNLS samples (5% for SNe la at z > 0.7), but 
still within the uncertainties. The decrease in fiducial 
rise-time with the 2-s respect to the l-s fit comes par- 
tially from the compensation for the higher rise-stretches 
measured. 

The 2-s technique introduces new parameters in the 
fit, and we want to assess if an improvement is sim- 



terion. BIC 



fln(N), a tool that gives weight 



to the number of para meters of the fit, /, as well as 
the total data used, N (Schwarz [1978f ). By summing all 
the individual core LC fits (— 10 < r < 35 days) of all 
SNe, we find that the l-s fit gives a smaller BICvalue: 
BIC = 17405, compared to the 2-s value: BIC = 19721 
(the average \ 2 P er degree of freedom for the individual 
2-s fits is slightly better, (xt) — 1-44, than for the l-s 

technique: (x 2 ) = 1-48), i. e. the 2-s technique did not 
lead to a better fit once the increase in parameters is 
taken into account. If we only consider the rise portion 
of the lightcurve (for which we use the same quadratic 
fit for both techniques) , we obtain a lower chi-square per 
degree of freedom for the l-s approach, xt = 1-60, than 
for the 2-s fit: \l = 1-64. The SNLS systematic uncer- 
tainties from different LC fits, 2-s and a newer SiFTO 
version different from Gil, are shown in Table [3j 

Similar results are found for the SDSS and low- 2; sam- 
ples. The core LC fits are better for the l-s fits whereas 
the early rise-time fit is better for the 2-s approach, ac- 
cording to the BIC criterion. The average rise-stretch 
is larger than the fall-stretch (by 0.09 for the SDSS and 
0.11 for the \ow-z sample) partly because of a mean shift 
in the peak date (of 0.45 and 0.40 days for the SDSS 
and low-z, respectively). The mean errors in peak time 
are somewhat lower (0.69 for the SDSS and 0.21 for the 
nearby sample). The final rise-times obtained are com- 
parable for both techniques, although the low-z 2-s result 
is ~ 1.5cr lower than its l-s counterpart. 

In the present study, we do not find a significant change 
in the rise-times of different samples for the two parame- 
terizations considered: l-s and 2-s, except for the nearby 
set for which a ~ 1.5er difference is found. Our 2-s tech- 
nique that is based on an adapted two-piece SiFTO fitter 
provides more uncertain peak date estimates. Requiring 
accurate maximum epochs affects the number of SNe we 
can use in the fit for the 2-s technique (see section [3]). 
We find that using only SNe with At$ < 1.2 days is good 
enough for the fit but the number of objects used is less 
than for the l-s technique, even though we use a tighter 
criterium for that fitter of At$ < 0.7 days. 

Although the statistical analysis of the fiducial rise- 
time remains unchanged with the parameterization used, 
Figure [4] shows that individual rise- and fall-stretches in 
the 2-s technique do not vary equally for different SNe. 
The linear fits show that lower-s SNe la tend to have 
larger rise-stretches and shorter fall-stretches for their 
respective l-s (a higher difference between the two) as 
compared to higher-s SNe la. The slopes of rise- and 
fall-stretch are different by more than 3cr. 



4.3. Comparison of different templates 
The temp l ate used in C06 descended 



from 

iNugent et al.l (|2002f ) and the early portion was based 
on a limited number of SNe la. Using this template 
to calculate the rise-time, we obtain 19.33 days, in 
agreement with the result obtained in C06 of 19.10 days. 

The fiducial stretch-corrected SN la rise-time is depen- 
dent on the SED template used (see Figure [6]), in par- 
ticular on the rather arbitrary definition of a s = 1 SN. 
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TABLE 2 

Stretch-corrected rise-time fit results for the single- and two- stretch lightcurve fits (with rt = — 10 days) 



Sample 


N SN {l-s) 


t r (l-s) (days)* 


iV S iv(2-s) 


t r (2-s) (days)* 


SNLS z < 0.7 


134 


17 02+ 19 


87 


16.92lO;f 9 


s < 0.8 


12 


18.03+°/ " 


8 


17.89+ij* 


0.8 < s < 1.0 
s > 1.0 


55 
67 


17.44«if 3 
16.87«;« 


31 

48 


I7.56±°f 7 
16 62+ 44 


SNLS z > 0.7 


102 


17 11+ - 24 


15 


16 11+ 068 


Nearby 


16 


17.39±££ 


15 


16.58io'27 


SDSS-II 


64 


17.44l°;i 


36 


17.03± ';f 


0.8 < s < 1.0 


23 


1 7 78+0-34 


12 


i7-5iiS:S 


s > 1.0 


41 


17 07+0.44 

-0.44 


24 





Statistical MC errors only 



TABLE 3 

Stretch-corrected rise-time systematic uncertainties for the SNLS stretch samples: low-s (s < 0.8), mid-s (0.8 < s < 1.0) 

AND HIGH-S (s > 0.8) AND DIFFERENCES BETWEEN THEM 



Systematic 




Rise-time t r (days) 




Difference At r (days) 




all 


low-s 
(s < 0.8) 


mid-s 
(0.8 < s < 1.0) 


high-s 
(s > 1.0) 


low-s - mids 


mid-s - high-s 


STANDARD 


17 02+ ' 19 


18.03t°J 7 


17 44+0.45 
11 '^-0.43 


16.87i°;« 


o.59±?;f 4 


57+0- 47 
U -°'-0.50 


LC-FIT 














1-s new SiFTO 


17.14 


18.46 


17.99 


16.57 


0.47 


1.42 


2-s 


16.92 


17.89 


17.56 


16.62 


0.97 


0.94 


TEMPLATE 














NEW (conley09f) 


16.94 


18.33 


17.70 


16.59 


0.63 


1.11 


C06 (nugent05)* 


18.77 


20.13 


19.68 


18.21 


0.45 


1.47 


SALT2 


18.16 


19.13 


18.77 


17.57 


0.36 


1.20 


TRANSITION EPOCH 














n = -8d 


17.22 


18.68 


17.56 


16.95 


1.12 


0.61 


PHOTOMETRY 














Fr phot, b = 


17.21 


17.96 


17.40 


16.96 


1.02 


0.44 


Ca phot, free b 


17.18 


18.07 


17.50 


16.86 


1.03 


0.64 


POWER LAW 














n = 1.9 


16.83 


17.87 


17.26 


16.68 


0.61 


0.58 


_tot 
CT syst 


+ 1.78 
— 0.23 


+ 2.24 
— 0.22 


+ 2.31 
— 0.18 


+ 1.35 
— 0.45 


+ 0.69 
— 0.27 


+ 1.24 
— 0.13 



This template is based on an old limited SN sample and is not included in the final systematic error budget. 



of the SNLS, like the one used in Gil, or a more recent 
version trained with a larger sample (used for SiFTO in 
IGuv et al.ll201(l . give consistent rise-times after correc- 
tion for the different s = 1 definition. Table [3] shows the 
rise-times for different templates. 

4.4. Comparison with SALT2 

We in vestigate here th e use of a different LC fitter, 
SALT2 (|Guv et all 120071) . This fitter also makes use of 
different templates from the ones presented in the previ- 
ous section, so that we are actually probing two system- 
atics simultaneously: template and LC fitter. 

In order to do this, we fit all SNe la with SALT2 
in the range r > — 10 days. A direct comparison of 
their SiFTO stretch and X\ parameter from SALT2 is 
shown in Figure [5] through a linear fit. We calculate 
restframe, stretch-corrected epochs using the redshift 
and a corresponding stretch (for each XI parameter): 
t = (t — t )/s(Xi)(l + z). Every data point matching 



This is especially important when comparing rise-times 
from different sources, where various templates can lead 
to different rise-times. A s = 1 SN is defined here as the 
mean stretch of the training sample for the template we 
use, which is the same as in Gil. If we correct rise-time 
of the C06 template for the different s = 1 definitions, by 
comparing the fits directly (see Figure!?]), i< e - a s = 1 SN 
for the Gil template corresponds to a s = 1.03 for the 
C06 template, we then obtain 18.77 days. This averaged 
correction partially accounts for the early wide behavior 
of the C06 template (see Figure [2]). The simple change 
of s = 1 definition is not sufficient to fully compensate 
the effect of an incorrect early template, so that a longer 
stretch-corrected rise-time is still obtained. We get sim- 
ilar long rise-times with the C06 template applied to the 
2-s technique, indicating that this parameterization is 
not independent of the template used. 
Newer templates trained with the well-sampled data 
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0.15 



0.8 1.0 
1 s - stretch 

Fig. 4. — Comparison of rise- (light blue triangles) and fall-stretch 
(dark blue inverted triangles) for SNLS SNe la calculated with 
the two-stretch technique as a function of the single-stretch. The 
crosses are the mean of the rise- and fall-stretches. The lines are 
linear fits: solid (red) for the mean of rise- and fall-stretches (slope: 
0.91 ±0.08), dotted (light blue) for rise-stretch (slope: 0.72 ±0.07) 
and dashed (dark blue) for the fall-stretch (slope: 1.09 ± 0.04). 




n L_n , n 

5 10 

to(S)-to(S„S,) 



-3 -2 

A t (s)-A t„(s„S,) 



Fig. 5. — Difference of peak date (upper) and peak date error 
(lower) calculated with the 1-s and the 2-s techniques for the SNLS 
sample 





Fig. 6. — Different SN la light-curve templates (inset of early 
portion on the right): C06 (solid black), Gil (dotted blue), NEW 
(dashed red) and SALT2 (dot-dashed green). 




Fig. 7. — Stretch comparison between lightcurve fits of the Gil 
template and other templates: C06 (red triangles), and a more 
recent one NEW (purple squares). The red dot-dashed and purple 
dashed lines are the averaged median in 0.05 stretch bins. 

restframe B-band is Jf-corrrected by t aking the SEP for 
each restframe epoch based on Eq. 1 of lGuv et al.l (|2007[ ) 
with Xq, X\ and the color obtained from the SALT2 fit, 
and integrating it in the appropriate filter at the proper 
redshift. We then fit the aggregate LC points of the early 
rise (r < — 10 days) to the quadratic rise model. 




Fig. 8.— Comparison of SiFTO stretch and SALT2 X\ for the 
SNLS SN la sample. The solid line represents the fit: st = axXi+b 
with a = 0.086 ± 0.026 and b = 0.997 ± 0.022. 

Similarly, we take into account errors with a MC that 
shifts the fit parameters (z p hot if necessary, and r max , 
Xq, X\ and color) according to their covariance matrix, 
and re-calculate a stretch from this, as well as a new K- 
correction for each point of each simulation. The median 
of all iteriations, the final fiducial rise-time we obtain for 
the normal- s population, is 18.16 ± 0.44 days. This is 
larger by ~ 1.1 days than the standard SiFTO 1-s and 
2-s rise-times. Although the SALT2 template is different 
than the templates we use (see Figure 0, it is not as 
extreme as the C06 template. Clearly, the use of the 
LC fitter also affects the measured rise-time and is an 
important systematic. The differences between SN la 
populations are still consistent as shown in Table [31 

4.5. Transition Epoch 

The transition epoch is the minimum epoch used in the 
core LC fit and the maximum epoch allowed in the early 
rise-time fit. Its value can affect the quality of the fit 
but also the resulting rise-time (and even the quadratic 
behavior of the rise). A higher value, like r t = —6 days, 
increases the amount of data for the early-time fit, but 



The Rise-Time of Subluminous SNe la 



9 



the rise may no longer follow a simple quadratic relation. 
Figure[9]shows that rise-time fits with different transition 
times and a fixed power fit (quadratic n = 2) for the 
SNLS do not change much up to a maximum value of 
r t ~ —8 days, where the rise-time starts to diverge. We 
use here the standard r t = — 10 days used in other studies 
but add a systematic uncertainty to the rise-time based 
on the differences at r t = —8 days (see Table [3]). 



. 1 s-f it 
2s-fit 



S 17.5- 
1 



Transition Time t, (days) 



Fig. 9. — Stretch-corrected rise-time as a function of the transi- 
tion epoch, the maximum date allowed in the single- (black circles) 
and two-stretch (red squares) rise-tim e fit s. The rise-time errors 
are the MC errors calculated as in § 14.11 The larger points at 
rt = —10 days are the quoted fiducial rise-times of this study. 



4.6. Zero scale b 

Wrong photometric subtraction affects fluxes that are 
close to zero and strongly influences the final stretch- 
corrected rise-time calculation. To study this systematic 
error, we add a nuisance parameter, b, to the quadratic 
form (Eq. [5]) that will absorb any general photometric 
error. For a good photometry, we should obtain b values 
close to zero. In our study, we get b = (— 0.48l ; 09 ) x 
10~ 3 , an off-zero value reflecting the large dispersion of 
the fluxes near zero (median of 3.1 x 1CT 4 with 0.03 
RMS) but small enough to represent a negligible effect 
on the rise-time. Consistent rise-times are obtained with 
free/fixed b values for the photometry processed via the 
French and Canadian reconstructions, as seen in Table 

El 

The import ance of this parameter can be observed in 
Figure 10(a)| in which different b values in the range 
—0.05 < b < 0.05, or a 5% maximum shift of the normal- 
ized flux at peak, result in rise-times shifted by several 
days. We emphasize therefore the importance of proper 
photometry, well-sampled to dates even before explosion, 
to correctly calculate 6 and not affect the rise-time. 

4.7. Different rise power laws 

Other parameters that could influence the fiducial rise- 
time behavior are a and the rise power n in Eq.[5J We find 
a = (7.40lg;^) x 10~ 3 for the normal-s SNLS population 
(arbitrary flux units) and an assumed n = 2 quadratic 
form. If we relax the quadratic rise assumption and allow 
the power of the rise-time model to be free (free a = 
(11.62 ± 0.85) x 10" 3 and fixed b = 0), we obtain n = 
1. 92^0.37 (stat). Rise-times fitted with this new power 



law are lower for all stretch populations but consistent 
with thei r quad ratic equivalents. 

Figure 10(b) shows the relation between the assumed 
power-law, n, and the measured stretch-corrected rise- 
time, t r for the normal SN la sample. A higher n makes 
the LC rise faster. In order to compensate for this, 
the duration of the rise to maximum of the fit becomes 
longer. This effect has a high degree of degeneracy and 
correlation between the two fit parameters, t r and n. A 
very similar trend is observed for the different stretch 
samples. Furthermore, if the quadratic form is relaxed 
for all populations, we obtain increasing power laws and 
rise time for higher- s SNe (Table 01 . This variation 
shows again the clear difference in the rise-time region 
for the sub-samples of SNe la. 



TABLE 4 

Stretch-corrected rise-time and power-exponents fit 
results for the different stretch samples in the snls 



Sample 


t r (days) 


* 

n 


SNLS 
High-s 
Mid-s 
Low-s 


16.85±°;£ 

16.49l°;|l 
16.141*;°° 


1 O9+0-31 

2 66+ 44 

l.48l°; 3 3 * 
0.78l°;£ 



Statistical MC errors 

Although the systematic uncertainty from this varia- 
tion is also highly correlated with the transition time, we 
still include it in the total systematic error budget. 

5. RESULTS AND DISCUSSION 

The final fiducial rise-time, t r , for all SNe la in the 
SNLS calculated via the single-stretch method is t r = 
17.02lo 26(statistical)l u ;23(systematic). The statistical 
uncertainty represents the lex distribution of the MC sim- 
ulations and the systematic uncertainty comes from the 
different effects considered in the previous section (added 
in quadrature) and shown in Table [3l This agrees with 
the value found via the two-stretch technique. 

We find that the fiducial stretch-corrected rise-times 
are consistent for fast and slow decliners, independently 
of different fit characteristics such as the template or 
lightcurve fitter. Mid-stretch SNe la (0.8 < s < 1.0) 
have a fiducial rise-time At r = 0.57lo^o( s t a t)lo;i3(syst) 
days longer than their high-stretch counterpart (s > 1.0), 
barely a ~ la result. With even lower significance, low- 
stretch SNe la show a difference with respect to the high- 

(stat)lJ;g(syHt) 



x - iu — 1.39^ 



stretch population of At r 
days. 

Although we do not measure individual rise-times, an 
estimate based on the stretch of each object can be ob- 
tained from the fiducial rise-time. For the single-stretch 
case we do r r .; = Sj x t r (s), where the subscript i de- 
notes each individual SN and t r (s) is the fiducial rise- 
time. Individual rise-times for the low-s population lie 
in the range of 10-14 days. For the normal-s population, 
they lie in the range of 14-21 days. 

5.1. Comparison to previous studies 

As mentioned in § 14.31 the different definition of a fidu- 
cial s = 1 lightcurve across different templates makes the 
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Fig. 10. — Rise-time t r as a function of zero scale b (upper) and 
power exponent n (lower) of Equation [2] as obtained in the fits 
for 1-s (black circles) and 2-s (red squares) tech niques (at r< = —8 
days). Errors indicate MC uncertainties (§ 14.11 1, Empty points are 
the best fits with b and n free, respectively (a was kept free for all 
fits). 



comparison of stretch-corrected rise-times from several 
studies quite difficult. 

Having that in mi nd, recent studies (jStrovinkl 120071 : 
lHavden et al.l l20~10bD . including the present one, agree 
that the average rise-time of SNe la is shorter than pre- 
viously thought. With more data, better-sampled at 
early times, the new templates, as well as LC-fitting tech- 
niques, provide evidence for a ~l-2 days shift in the av- 
era ge rise-time of SNe la with respect to previous stud- 
ies flRiilLjEII] EH1 iConlev et all [20061. although the 
precise value is heavily dependent on the definition of a 
fiducial s — 1 LC template. 

Although we use the same technique as C06, the tem- 
plate they used was trained with a handful of SNe la and 
consequently had a rise-time behavior to maximum quite 
different from the current templates. Their rise-time is 
longer than the one we find here, even after correcting 
for our s = 1 definition. 

We also analyze the nearby and SDSS-II datasets with 
our same template and LC fitter. H10 used a much larger 
SDSS-II sa mple (361 SNe) than us ed here (the 146 SDSS- 
II SNe from lHoltzman et al.ll2008h . Our SDSS-II fiducial 
rise-time agrees well with theirs. They argue that the 
two-stretch technique presents a much better parameter- 
ization than a single parameter approach. We do not 



agree with those findings for our aggregate rise-time fit 
approach, neither for the SNLS nor the SDSS and low- z 
data. This could be due to the use of a different fidu- 
cial LC template and fitter that differ in the rise-time re- 
gion from t he ones they us e (Gil and modified SiFTO vs 
MLCS2k2. |jha et al.ll2007D . in particular the stretch they 
calculate is not based on multiple bands. The Bayesian 
Information Criterion for the core LC fits suggest that 
the inclusion of an additional parameter to model SN la 
lightcurves does not provide a valuable improvement of 
the fit. The quadratic fit of the early rise-time data is 
not strongly affected by the parameterization either. In- 
dividual rise-times, on the other hand, could benefit from 
such an independent parameterization, as shown in Fig- 
ure m 

Based on our calculations of the fiducial rise-time for 
different samples across the redshift range < z < 1, 
we do not find evidence for a significant evolution in the 
early portion of SN la lightcurves. Dividing the SNLS 
sample into three populations according to their stretch 
values, we measure the average rise-times for each sub- 
group and find them to be consistent. Although at low 
confidence, the data also agree with a trend of faster fidu- 
cial rise-times for brighter and slowly declining SNe la 
after the primary stretch correction, for the SNLS and 
SDSS-II. H10 obtained a similar trend for the individual 
rise-times of SDSS-II objects. 
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Fig. 11. — Distribution 56 Ni masses (solar units) for the normal- s 
(solid histograms) and low-s (filled histograms) SN la SNLS sam- 
ples based on 1-s (black) and 2-s (red) fits. Vertical lines represent 
the mean for each case. 



5.2. 5e Ni masses 

In order to estimate the 56 Ni masses, we need to cal- 
culate bolometric luminosities. We integrate the model 
SEDs for each SN across all observed bands. This omits 
the near- infrared (NIR), which accounts for a substantial 
fraction of the bolomet ric luminosity , so w e must apply 
a correction factor. In IHowell et al.l (|2009Q the IR per- 
centage for normal SNe la is estimated via integration of 
proper SEDs corrected for each SN in the region outside 
the observed spectrum, under some assumptions. For 
low-s SNe la, SEDs are not as fully sampled. We instead 
use a very broad approximation based on the IR frac- 
tion found for SN1999by and SN2004eo (IHoflich et al.1 
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driver, when studying early-time effects of SNe la. 




Fig. 12.— 56 Ni mass as a function of individual rise-time for 
the single- (upper) and two-stretch (lower) techniques. The low- 
s (orange circles), mid-s (blue triangles) and high-s (light blue 
squares) SNe la are indicated. 



l200l iGarnavich et~aI|[2CJ0l iTaubenberger et al.ll2008h : 
~ 15 — 25%. The fraction is at least 5% higher than for 
normal SNe la. We use a standard 20 ± 10% value (and 
10 ± 5% for the normal-s) to calculate via Equations [3] 
and |4] the 56 Ni mass. 

In order to account for the errors in the LC parameters, 
we carry out this computa tion for every realization of 
the LC data in the MC of § 14. 1[ and use the median Lb 
of each SN to calculate the final 56 Ni masses, as shown 
in Figures [11] and [T^l We also include the error on the 
conversion factor 7, which can also slightly vary with 56 Ni 
mass. The 56 Ni masses lie between 0.05-0.35 Mq for the 
low-s and 0.1-0.8 M Q for the normal-s populations. Our 
estimates are lower than typical values in the literature 
because of the lower rise-times. They are similar for both 
parameterizations of the lightcurve (see Figure [T3")) . 

Figure [12] shows 56 Ni masses as a function of the indi- 
vidual rise-times for all SN la. The relation of Ni mass 
and rise-time is expected from the use of Arnett's rule. 
The individual rise-times are obtained from the fiducial 
rise-time and the individual stretches, so that 56 Ni mass 
also relates with stretch. However, one can observe that, 
as rise-time is more independent of stretch for the 2-s 
technique, there is a more diverse 56 Ni mass range for 
each stretch population. It may be beneficial to divide 
the LC into separate regions, rise and fall, to disentan- 
gle the rise portion from the primary width-luminosity 
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Fig. 13. — Comparison 56 Ni mass estimates as calculated in § 15.21 
for the single- and two-stretch techniques. Orange circles are the 
low-s objects. 



5.3. Secondary lightcurve parameters 

The need for secondary parameters to describe SNe la 
at early times, beyond the traditional core LC parame- 
ter like stretch or A77J15, would have important conse- 
quences for the understanding of SNe la. To the preci- 
sion achieved in this study, these secondary effects in the 
early rise-time region do not seem necessary to describe 
SN la LCs. We find a low-significance trend towards 
longer rise-times in fast decliners after stretch correc- 
tion. If confirmed with more data, this would mean that 
rise-times do not correlate to second-order with stretch 
and 56 Ni mass. This would suggest that beyond the ra- 
dioactive decay of nickel which only sets in later in the 
LC, some other parameter in the initial explosion in the 
outer ejecta could play a role. 

Besides a primary parameter responsible for the 
luminosity-width relation of SNe la, like 56 Ni mass (from 
a varying density at the transition of subsonic deflagra- 
tion to su personic detonation in the delayed-detonation 
scenario). IHoflich et al.l (|2010D predict that there could 
be other variables like the C/O ratio of the WD progen- 
itor that would affect the early part of the lightcurve. 
Higher C/O ratios result in a smaller binding energy 
at time of the explosion of the WD. In turn, this leads 
to higher expansion velocities meaning a faster photo- 
sphere recession and a shorter rise to maximum light. 
Lower main-sequence masses (during He-core and He- 
shell burning for example) can lead to higher C/O ratios. 
This would imply that lower-stretch SNe la with longer 
rise-times come from higher masses. The fact that low- 
stretch SNe la are primordially found in old environments 
sugges t that their progenitors are l ong-lived and of low- 
mass fGonzale z-Gaitan et al.l [201 If) . in contradiction to 
these findings. 

On the other hand, low-stretch SNe having longer 
rise-times is remi niscent of the behavior predicted by 
iFink et all (|20100 for the double-detonation of sub- 
Chandrasekhar models. The initial He-shell detonation 
on the WD surface can lead to a detonation of the CO 
core of the WD. Varying WD masses (0.8-1.4M Q ), and 
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with it their central core and shell densities, give rise to 
a variety of 56 Ni masses and burning products power- 
ing the lightcurve. Particularly they find that brighter 
events have faster rise-times. This is because brighter 
explosions produce more core 56 Ni that expand faster 
with less opacity in the outer shells from iron-group- 
elements created in the He-shell detonation. Their ef- 
fect is anti-corrected with stretch, which would indicate 
a primary correction rather than a secondary one. In 
this picture, the lightcurve and the rise-time are directly 
related to the initial mass of the WD. Lower-mass WDs 
(of ~ 0.8Mq) have lower He-shell and CO core densities 
leading to a distribution of burning products that give 
rise to slower rise-times for low-stretch SNe la. Lower 
mass WDs producing fast decliners would agree with the 
aforementioned observations of older host environments. 

Another intere sting feature of the predictions by 
iFink et abl (|2010D is that the contribution of other ra- 
dioactive elements like 52 Fe and 48 Cr of short decay times 
can make the rise-time longer. This is particularly true 
for subluminous SNe la, which in their models produce 
as much as 20% of these radioactive nuclei, particularly 
44 Ti, which does not contribute to the early lightcurve, 
but is known to be observed in the spectra as well. These 
yields are a direct consequence of the initial He-shell det- 
onation. 

6. CONCLUSIONS 

We have constructed aggregate restframe i?-band LCs 
of a sample of spectroscopic and photometrically identi- 
fied SNe la at z < 0.7 from the SNLS. We obtain fiducial 
rise-times that are consistent with measurements for the 
low-z and SDSS-II samples calculated here with the same 
techniques. We extend our SNLS sample and show that 
there is no significant evolution in the early portion of 
SN la lightcurves up to z = 1.0. 

We find that the use of a two-stretch parameteriza- 
tion for the lightcurve fits of SNe la is equivalent to the 
single-stretch parameterization for the rise-time study of 
aggregate lightcurves. We confirm the consistency for 
the SDSS-II and SNLS samples while the nearby set 
has a rise-time ~ 1.5a lower for the 2-s technique and 
discrepant with all other calculated rise-times. We in- 
sist in the importance of proper lightcurve fitters with 
parameters that can take into account earlier parts of 
the lightcurve (stretch vs Amis) and that include multi- 
band information, and of templates trained with suffi- 
cient SNe with early time data. We also find indications 
that the 2-s approach might be useful for detai l ed stud - 
ies of individual rise-times, as in lHavden et al.l ((2010b); 
Gancshaling am et al.l (|2011l ). 

We have compared the early region fits for different 
stretch population of the SNLS. We find that the stretch- 
corrected rise-times are consistent within uncertainties 
for slow and fast decliners. We extend this result to the 
lowest stretch objects. These findings are independent of 
systematics such as the maximum epoch allowed in the 
early time fits, as well as on the two LC parameteriza- 
tions assumed. The present study averages over samples 
which can have obvious difficulties like the low statistics 
of the low-stretch population. At a low significance of 
~ lcr, our results also agree with a trend for longer rise- 
times in fast declining SNe la. If these indications are 
confirmed, then SN la explosions differ beyond a single 



parameter orchestrating the lightcurve of SNe la. Sec- 
ondary parameters in the early rise region would be nec- 
essary to fully model the lightcurves and could provide 
useful insights into SN la explosions. Such an indication 
is also found for the SDSS-II, as in H10. 

However, to the precision attained in this study, SNe la 
remain extremely homogeneous all the way to very early 
times. One single lightcurve parameter, that possibly 
translates to a unique physical variable, like the tran- 
sition density of deflagration to detonation in the de- 
layed detonation model, can account not only for the 
core lightcurves of SNe la, but largely for the rise-time 
as well. The similarity of SN la lightcurves (after stretch- 
correcting) argue for a real common origin, with sublu- 
minous SNe la at the extreme end of the population. 

The fiducial rise-time we find translates into 56 Ni 
masses of O.O5-O.9OM for all SNe la (and 0.05-0.35M o 
for low-s SNe la). Although our assumed IR contri- 
bution to the bolometric luminosity is uncertain, ex- 
tremely low amounts of radioactive 56 Ni seem enough 
to power a faint SN la, l es s tha n that the amount 
found by IFilippenko et al.l ((1992D ; ILeibundgut et al.l 
(1993). Our SNe with the lowest radioactive nickel yield 
even rival wit h peculiar extrem e obje c ts such as the 
ones fo und by iMcClelland et all (j2010fi : iKasliwal et~aH 
((201(1 ; [Poznanski et all (|2010l) . and pose a challenge to 
explosion models. 

Such weak explosions seem unlikely to be the outcome 
of a runaway of a WD near the Chandrasekhar mass. 
Sub-Chandrasekhar models have been indeed a natu- 
ral way of explaining sublumi nous SNe la (e. g. iLivnd 
[19901: IWooslev fc Weaven[l99l. an d possibly all SNe la 
(IFink et all 120101: iSim et al.l l2010t Ivan Kerkwiik et al.l 
120101). Located in passive environments that are charac- 
teristic of old populations, low-mass WDs seem a promis- 
ing progenitor of subluminous SNe la. Wit h recent rejec- 
tions of a high fraction of r ed giant donors (jHavden et al.l 
I2010bt iBianco et al.ll20rll ). a double degenerate merger 
leading to the sub-Chandrasekhar explosion of a low- 
mass WD with very low 56 Ni yield appears increasingly 
plausible. If SNe la share a common explosion and 
progenitor mechanism, a varying progenitor WD mass 
could possibly explain the different SN la radioactive 
yield and luminosity, and their environment. More and 
better simulations n e ed to test the validity of this pic- 
ture. Pakmor et al. (2010), for example, modelled the 
merger of two ~ 0.9M Q WDs and obtained a sublumi- 
nous SN la explosion with ~ 0.1M Q of synthetized 56 Ni, 
in agreement with the bulk of our low-s objects but in 
disagreement with the expected progenitor masses from 
the environments. More theoretical development need to 
complement the abundant current and coming observa- 
tions of SNe la. 
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